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ABSTRACT 


i 


The effects of vibrational and chemical nonequilibrium on 
turbine performance were investigated separately. The vibra- 
tional model was taken as a pure nitrogen expansion, and the 
chemical model was taken as the combustion products of a 
stoichiometric mixture of kerosine and air. The loss of per- 
formance of fully frozen flow with respect to the equilibrium 
flow for each model was determined. The extent to which non- 
equilibrium will occur was investigated within limited ranges 
of pressure and temperature. Vibrational nonequilibrium can 
result in losses up to seven percent with respect to equilib- 
rium flow. Vibrational freezing will be virtually complete a 
short distance after the throat of the first stator assembly. 
The losses due to chemical nonequilibrium are insignificant 


compared to the vibrational nonequilibrium losses. 
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1. INTRODUCTION 


The aircraft of tomorrow must be capable of operating at 
high Mach number, have large amounts of excess thrust, be 
capable of operating at high altitude, and yet be very eco- 
nomical. As the requirement for these aircraft becomes more 
prevalent, designers of both airframes and power plants must 
strive to match their respective ideas to produce a viable 
design. The airframe specialist is concerned with drag and 
high skin-temperature problems which he hopes to resolve by 
reduction of frontal area and by lighter gross weights. The 
power plant specialist investigates higher turbine inlet tem- 
peratures because of the increased efficiency of the cycle 
with high peak temperatures. The increase in flight Mach 
number induces a rise in the temperature of the cooling air 
available at the turbine. This higher temperature reduces the 
cooling potential of the air. 

The next generations of the turbojet engine will utilize 
nearly double the present turbine inlet temperature by burn- 
ing stoichiometrically. As the turbine inlet temperature is 
increased, the effects of dissociation and vibrational excita- 
tion play a significant role in the energy balance. The jump 
in turbine temperature for stoichiometric combustion will 
produce almost triple the degree of vibrational excitation, a 
large increase in local velocities for the same local Mach 
number within the engine, and a substantial rise in the degree 


of dissociation of the products of combustion. 
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The phenomenon of vibrational energy freezing occurs due 
to the relatively large number of molecular collisions re- 
quired to return the vibrational mode of internal energy to 
equilibrium. Freezing of chemical composition occurs when 
the gas cannot react fast enough to maintain equilibrium with 
a rapidly changing temperature and pressure. Both of these 
phenomena are familiar to the rocket designer. The relatively 
low temperatures and high pressures characteristic to the 
turbojet cycle in the — have not caused these nonequilib- 
rium effects to any significant degree. 

The presence of nonequilibrium in the form of frozen 
vibrational energy and/or frozen chemical composition is very 
probable in the turbine section of an advanced high tempera- 
ture turbojet. The rapid removal of work from the gas in the 
turbine reduces both the temperature and pressure at a rapid 
rate. The vibrational temperature is the independent variable 
in the vibrational energy distribution, and there is a cor- 
responding temperature for both rotation and translation. At 
equilibrium, the rotation, translation, and vibration tempera- 
tures are the same. Comparison of the vibrational, rotational, 
and translational temperatures in the turbine would show a more 
rapid decrease in the latter two relative to the rate of 
vibrational temperature decrease, primarily due to the number 
of collisions required to maintain equilibrium. The magnitude 
of the relaxation time relative to the short transit time in 
the turbine suggests that a significant degree of vibrational 


freezing may result. The effect of lower pressure would be to 


Tigh 





reduce the collision frequency, thereby increasing the relaxa- 
tion times. 

Evaluation of the loss of available energy due to 
vibrational freezing becomes rather complicated in multi- 
component gases and involves complex-numerical computer pro- 
grams. Stollery and Park [17] and Wilson, et al. [22] give 
examples of these programs. Exact calculations of recombina- 
tion and dissociation rates of a mixture of gases such as 
combustion products also involve complex computer programs. 
The complexity of the entire problem is a result of the num- 
ber of nonlinear equations involved and the lack of an exact 
solution to the flow problem in rotating machines. 

This thesis attempts to apply the nonequilibrium analysis 
of vibrational and chemical freezing to the turbine section 
of a turbojet cycle. The presence and degree of each of these 
problems are investigated, and the effects of changing the 


thermodynamic parameters of the cycle are discussed. 
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ii.) ANAS £5 


High velocity flows frequently create rapid changes in 
the flow variables: temperature, pressure, and density. The 
internal modes of freedom of the gas molecules do not always 
remain in equilibrium or return to equilibrium fast enough to 
keep pace with the rapidly changing macroscopic environment. 
The return to equilibrium of these internal modes is governed 
by internal kinetic equations. The existence of chemical re- 
actions in the flow introduces the chemical kinetic equations. 
There are characteristic times associated with each of these 
processes: relaxation time with the internal kinetic equa- 
tions and reaction times with the chemical kinetic equations. 

Maintenance of equilibrium of chemical reaction and of 
internal degrees of freedom requires a number of molecular 
collisions. It is obvious that true equilibrium flow cannot 
exist. The average time for a unit mass of gas to transit a 
turbine is a prime factor in determining whether or not 
freezing occurs. The conversion of available energy to useful 
work by the turbine is accomplished in a small distance, and 
present turbine technology involves Mach numbers in the 
transonic range. These facts, combined with high local tem- 
peratures, produce a local fluid velocity relative to the 
blade of the order of two to three thousand feet per second. 
The short physical length of the turbine results in an ex- 


tremely low transit time. If the characteristic times of 
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progression of the chemical and molecular phenomena are very 
small compared to the time intervals of interest in the flow, 
the flow ats be considered to be in ‘local equilibrium. Con- 
versely, 1f these characteristic times are much larger than 
the time intervals of interest in the flow, the phenomena 
May be considered frozen. 

Theoretical description of a nonequilibrium flow requires 
incorporation of the kinetic equations describing the processes 
into the usual gas-dynamic conservation equations. Analysis 
and solution of such systems of esac ave dilfti cule.) an 
all but the simplest problems, numerical methods are required. 
Bray and Appleton [5] and Anderson, et al. [2], have presented 
detailed analytical solutions to these problems. 

Extensive research into nonequilibrium flows was inspired 
by the advent of space exploration. The flow around re-entry 
vehicles at high altitudes involves complex nonequilibrium 
chemistry [6]. In rocket propulsion systems, the problems of 
large expansion ratios and low exit pressures tended to magnify 
the situation. Excellent literature surveys by Read [15], Rich 
and Treanor [16], and Taylor and Bitterman [18] provide a 
myriad of data on vibrational relaxation. Gas dynamic lasers 
employ intentional vibrational nonequilibrium to obtain popula- 
tion inversions [3]. 

The complexity of solution of the expansion of combustion 
products through a turbine should be readily apparent. The 
problem must therefore be simplified in order to obtain 
meaningful conclusions as to the effects of nonequilibrium 


Phenomena on the expansion. The interaction of vibrational 
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excitation and chemical reaction is not solved here, but some 
insight into the level of existence of each phenomenon is 
presented. 

The effects and mechanics of vibrational nonequilibrium 
flow are considered separately from the effects and mechanics 


of chemical nonequilibrium flow. 


A. VIBRATIONAL NONEQUILIBRIUM FLOW 

In order to investigate the effects of vibrational non- 
equilibrium in a turbine, a simplified model was utilized. 
The expansion of pure gaseous nitrogen through a turbine was 
considered instead of the more complex case of the actual 
combustion products. Nitrogen was a good model for a number 
of reasons. By uSing a pure gas, no chemical reaction, other 
than dissociation, was present. The degree of dissociation 
in the temperature range considered was negligible. The prob- 
lem of vibrational freezing could therefore be isolated, and 
the calculation of thermodynamic data for the flow was greatly 
Simplified. In addition, the mole fraction of nitrogen is 
large in the gas mixture. 

The following assumptions were expeditious: 

1. Complete (microscopic and macroscopic) equilib- 
rium existed at the exit of the combustion 
chamber. 

2. The gas behaved as an ideal, inviscid fluid. 

3. The turbine expansion was adiabatic. 

4. The flow was considered as quasi one-dimensional. 

5. The velocity of the gas at the combustion chamber 


outlet was considered negligible compared to the 
gas velocity after the turbine exit. 


ils, 





A turbine schematic (Figure 1) was defined. At station 
one, corresponding to burner outlet, the nitrogen was con- 
Sidered to be in complete equilibrium. Between stations one 
and two, the gas wasS expanded isentropically through the 
turbine. The work extracted by the turbine was assumed to be 
drawn strictly from the translational and rotational modes of 
internal energy. These modes were considered to be in equi- 
librium with each other during the entire process due to the 
relatively small number of collisions required for equilibrium. 
This model assumes the independence of the separate modes of 
internal energy. The validity of this procedure may not be 
universal in its application but was utilized here for sim- 
plicity of analysis. The vibrational energy was considered 
frozen in this expansion process. (See Figure 2.) 

Between station two and station three, the gas was allowed 
to return to complete equilibrium in a constant area duct. 
This was accomplished by the redistribution of excess vibra- 
tional energy to the rotational-transational gas through a 
process equivalent to Rayleigh heating. The rotational and 
translational modes were considered to be always in complete 
equilibrium. 

The extraction of work results in a decrease of total 
enthalpy. A hypothetical stagnation temperature of the 
rotational-translational gas at station two, te! was chosen 
to represent a given amount of work extracted by the turbine. 
The prime represents the state of the rotational-translational 


gas. The ratio of specific heats, y, was taken as a constant 
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and equal to that obtained by classical equipartition con- 
cepts for a diatomic gas with only rotational and transla- 
tional modes. This value is y = 1.4. The selection of a 
burner outlet stagnation presSure and temperature and a Mach 
number at station two, based on the chosen rotational- 

ies ational (RT) gas temperature, enabled the determination 
of the station two stagnation pressure through the use of the 


ideal gas isentropic relationship in equation (1). 
selene 
oe Neat s 
P! = Pp —— (1) 
Use of an isentropic procedure is consistent with the desire 
to isolate vibrational nonequilibrium. 
Downstream of station two, the gas was subjected to a 
Rayleigh heating process due to the relaxation of the vibra- 


tional mode. 


(2) 


Ny] Cc 
IES 


Utilizing the definition of stagnation enthalpy, 
' a t 
Q. + Hy H (3) 


where Q. represents the energy added to the RT gas from the 


vibrational mode. 


Equation (3) may be rewritten as: 


Q = Hi -~- HA = C oes: TY 2 (4) 
1g 3 2 Pom t. t. 


The ideal gas isentropic relationship was used to find static 


temperature, Toe which was needed in the iteration procedure. 


Hy) 





Ife 


2 
i ln ae (5) 
Neale ee 
oct 5 M, 


The determination of the Mach number at point three was then 
required. This value was determined by an iterative process 


Matching the Rayleigh heat added to the classic definition of 


enthalpy. The enthalpy was considered as 


T 
h(t) = J Ch aT (6) 


bt 
O 


Since the internal energy modes were considered separable, the 


Specific heat may be expressed as 


C = C + + C — + Cc (7) 


T T | 
h(T) = 5 eee J a, ht (8) 
ens Vv 


The energy added to the equilibrium RT gas in the Rayleigh 
heating model must be identically equal to the difference in 
the vibrational energy of the complete gas between stations 


One and three. This difference was taken thusly 


I 
2 
- 
‘?) 
tg 
< 
Qu 
ae. 
+ 
rH 
1 
@ 
rg 
fe 
Qu 
| 


h(T,) (9) 
Ak di 
Oo Oo 
Therefore 
ba 
h(T,) = Cp at + Cp (T)-T.) 
m V RT 
Oo 
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Similarly 


T 


3 
h(T,) = \ Cy aT + Cp (T,-T)) (10) 
T V RT 


O 


Taking the difference between equations (9) and (10) 


a 
aL 
h(T,) - h(T,) = J Cp dT + Cp (T,-T3) (11) 
1. V RT 
3 
Now define, Qne as the difference in the vibrational energy 
between stations two and three. 
‘aa 
Qn = PS aT (23) 
me 


Substituting equations (11) into equation (12) 


Q = h(T 


D = h(T,) wee 


p__ (T,-T3) (12) 


) 
- RT 


The iteration was satisfied when on of equation (4) was equal 
Ee) Qn of equation (12). The actual solution for a given set 
of input parameters was obtained by entering gas tables [11] 


with the chosen M.. With the already chosen value of Ty , the 
Z 


values of Tf and T* were obtained. A Mach number at station 


3' was assumed, and equation (4) was solved in the following 


Porm , 








t t 
nee iS 
Q. = Te Cp = (14) 
* * 
RT Te Te 


Utilizing information from gas tables [11] the temperature, 
Tae was also obtained and equation (12) was evaluated. This 


process was repeated until Q.. and Qn were within 0.5% of each 





other. Satisfaction of the iteration process allowed calcula- 
tion of stagnation temperature and pressure at station three. 

The investigation of performance lost due to vibrational 
freezing required a reference quantity. This reference quan- 
tity was taken as the turbine work produced by the expansion 
Of the gas with full equilibrium between the modes of vibra- 
tion, rotation, and translation. For both equilibrium and 
nonequilibrium cases, the turbine pressure ratios were assumed 
identical. 

Efficiency was defined as the ratio of actual work to 
ideal work where the actual work was considered as the work 
obtained from the vibrationally frozen model and the ideal 
work as that obtained from an equilibrium gas expanding to the 
same stagnation pressure. 


vibrationally frozen work 


Ny ™ equilibrium work is) 


H(T, ai H(T, ) 


i 3 
nN oa H(T. ) - H(T*® ) 
Vv H(T, H(TF ) (16) 
uf 3 
where Te is the ideal stagnation temperature at the stagnation 


3 


pressure of station three. 
The stagnation temperature for an ideal turbine evaluated 
at the stagnation pressure of station three was determined by 


isentropic relationships. By definition, 


| PP. ) = S(T* , P* ) aus 
Sq, a Si Se) ee. 


S( LeeLee = SCE) y= Reen (2) (18) 


20 


where S°(T) is the entropy of the gas at a given temperature 
and a pressure of one atmosphere. Using equation (18) in 
equation (17) 
So(T)) =. RQ = $°(T* ) - R &n P* (19) 
ty ty t. C3 


It was chosen that P*¥ =P 


. t and the values of T and P 


3 3 oT eT 
were assumed. Therefore, an equation for the entropy at the 


desired temperature was obtained. 


px 
a t. 
S°(T* ) = S°(T, ) + R fn | (20) 
eS ic Ee 
3) IL te 
1 
Solution of this equation yielded Tf which was utilized to 
3 


compute values of efficiency from equation (16). This effi- 
ciency represents a loss of performance resulting from non- 
equilibrium flow. | | 

A method to identify the significance of vibrational non- 
equilibrium had been established, and the question of actual 
presence of this phenomenon in a high temperature turbine 
remained. 

The complexity of the problem of vibrational relaxation 
precludes actual solution of the theoretical equations within 
the scope of this thesis, and experimental results of other 
researchers will be utilized in the ultimate discussion of 
the effects of this phenomenon on turbine flow. For complete 
understanding it is prudent, however, to develop the princi- 
ples involved. 

The nitrogen gas was considered to be a system of harmonic 
oScillators. This system of diatomic molecules was assumed to 


be capable of collisional interchange of energy with the 
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translational and rotational degrees of freedom of a heat 
bath at temperature T. The heat bath was provided by the 
translational and rotational degrees of freedom of the pure 
_gas. At equilibrium, the distribution of oscillators over 
the energy states was assumed to be given by the Boltzmann 
distribution. 

The majority of data on vibrational nonequilibrium, both 
theoretical and experimental, has been developed for the lag 
of vibrational equilibrium when vibrational degrees of free- 
dom are being excited. Rich and Treanor [16] point out that 
there 1s evidence that the relaxation time for de-exitation 
may be one-fifth or smaller relative to the exitation relaxa- 
tion time. The relaxation time is a measure of the rate of 
return to equilibrium after a disturbance. 

The equation for the time rate of change of the number of 
molecules in a particular energy state may be developed using 
Statistical thermodynamics and quantum mechanics. The rate 
of transition from one state to another is dependent upon the 
number of collisions per unit time between the oscillators 
and the molecules of the heat bath that possess sufficient 
energy to cause transition and on the probability that a 
Single collision will result in transition. 

For gas dynamic purposes the total vibrational energy is 
more useful than the molecular energy level population in 
each energy level. This total energy is a summation over all 
possible levels of the product of the energy in a particular 


level and the number of molecules in that level. Manipulation 


Ze 


of the rate of change of the molecular level population equa- 
tion produces an equation for the rate of change of total 


vibrational energy. Vincenti and Kruger [20] express this as 


dE, ES ~ EY 
at ~ T ae 


where EL 1s total vibrational energy, Ex is total vibrational 
energy at equilibrium, and tT is the relaxation time. 

This equation shows that the vibrational energy of the 
system will always move toward equilibrium and that the rate 
of this movement at any instant is linearly related to the 
degree of departure from equilibrium at that instant. Even 
though this equation assumed a Boltzmann distribution only at 
equilibrium and no restriction is placed on the amount of 
departure of vibrational energy from equilibrium, Vincenti and 
Kruger [20] limit the equation to small departures from equi- 
librium. The ieee to small departures is a result of 
the anharmonicity of actual molecules. 

The relaxation time, t, is a function of the temperature 
and pressure of the heat bath and of the characteristics of 
the vibrating species. For practical purposes it was assumed 
that equation (21) was valid for any number of excited mole- 
cules and could be applied to a time varying heat bath as long 
as EX and t were taken as functions of the instantaneous state 
of the heat bath. Under these assumptions, t should be tech- 
nically referenced as a local relaxation time where numerical 
integration could produce a solution for a particular time 


span. 
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For sufficiently low temperatures the original Landau and 


Teller Theory [20] may be approximated by: 


K 1/3 
ee S exp (x) (22) 


where K 1S a constant dependent upon the physical properties 
of the species, and C is an empirical constant dependent upon 
the type of heat bath molecules. 

In the present theoretical state, the evaluation of T must 
be accomplished by experimental measurements interpreted in 
the light of the foregoing theoretical results. The methods 
of ultrasonic absorption, shock tube interferometry, spectro- 
scopic data, and pressure measurement, have been used experi- 
mentally to produce the values of the constants for use in 
equation (22). These values are presented in Table II. 

Taylor and Bitterman [18] present numerous graphs of both 
t and the probability per collision of an exchange of energy 
between the translational and vibrational modes of a molecule. 

Anderson, et al. [2], presents correlations between re- 
laxation times, temperature and pressure for nitrogen and 
Various collision partners. These correlations were utilized 


extensively and are therefore presented here. 


ai 1 
105y9 (TP) yn, = 93 (fT -/°) - 4.61 (23) 
10g19 Py, -co, = Logg (TP) yon, (24) 
= -1/3 
Lod yg (TP) Ho = 15.4 (T ee i 22 (25) 
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Bray [4] suggests that once an appreciable deviation from 
equilibrium has occurred that return to equilibrium ina 
nozzle is unlikely. Therefore, the investigation was cen- 
tered on the stator nozzle block of the turbine. A stator 
passage with area change from throat to exit given by equation 
(26) and shown in Figure 3 was chosen, 

A/A* = 0.15x - 0.01x? (26) 
where A* was the area at the throat and x was the distance 
along the mid-span streamline measured from the throat. While 
the choice of this particular formula for change of passage 
area was arbitrary and not particularly critical, this choice 
is representative of modern turbines. 

Knowledge of the area change and selection of a represen- 
tative y enabled determination of Mach number, and thermo- 
dynamic data from gas tables as presented in Table III. 
Selection of a total pressure and temperature at the burner 
outlet permitted evaluation of vibrational relaxation times 
at the throat of the stator with area change given by equation 
(26) for nitrogen to nitrogen collisions using equations (22), 
(23), or the graphs of Taylor and Bitterman [18]. (See Table 
BY.) 

As stated previously, the calculated relaxation times had 
to be compared to a gas transit time. The gas transit (dwell) 
times were determined by dividing a characteristic length by 
the appropriate velocity. Under this concept, a total transit 
time was obtained by dividing the distance from point A to 


point B along the mid-span streamline by the mean velocity 
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between those two points. A local dwell time at a data point 
was obtained by dividing the distance between data points by 
the velocity at the upstream point. 

Numerical evaluation of equation (21) was performed 
utilizing calculated relaxation times and the method of Hurle, 
Russo, and Hall [9]. The method is illustrated in Appendix C. 
This evaluation illustrated the actual amount of vibrational 
energy as compared to the equilibrium amount. 

An aircraft gas turbine using hydrocarbon fuels will have 
Many molecular species in the combustion products. It is 
known that some of these species are more effective for vibra- 
tional de-excitation. Consequently, a composite relaxation 
time involving several molecules was calculated. The molecular 
Species in the combustion products were screened to ascertain 
the molecules most effective for vibrational de-excitation. 


This screening indicated HO and CO, should be included. A 


Z 
relaxation time was thus calculated for nitrogen with three 
different collision partners, No, CO,, HO, uSing equations 
(23), (24), and (25) respectively (Table V). The total rate 
of change of vibrational energy of the nitrogen was considered 
to be the sum of three distinct quantities corresponding to 
the change induced by each of the collision partners. 

dE dE. dE dE, 

ate otal aN -N . SAN, -H,0 iil 
Fach of the terms on the right-hand side of this equation was 


expressed in the form of equation (21) and substituted to form 


equation (28) 
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ay Ey Bo By ey 
total N,-N, N,-CO, N,-H,0 
The coefficients (A, B, C) were taken to be the mole fractions 


of the collision partner. Since the relaxation times with Ny 


and CO. are equal, equation (28) may be expressed as 


dE 1 - Te “HO 
ec ee 
V V 


a (29) 


where Xi 0 is the mole fraction of water in a gas mixture. 
2 


Comparing equation (29) with equation (21) 


i - *u,0 Ro 


L a 
saa cieeaaiada TE ean (30a) 


* composite "Ny Ny N.-H,O 


The extension of equation (30a) to systems with a larger num- 
ber of components is obvious. Note that EL, is the vibrational 
energy of nitrogen only. Values of t composite are in Table 
wr. 

Consider a large number of NX, where NZ denotes vibra- 
tionally excited N,- If the number of N%X¥ exceeds the equilib- 


2 


rium value, there will be a decrease in the number of N5 as 


time progresses. The N3 may be de-excited by collisions with 


HO iG CO. or N.- The collision frequency with H.,O, CO, , Or 


2 
N, depends upon the mole fraction. The probability that a 
given collision will result in de-excitation depends on the 
partners; i.e. on C and K, in equation (22). Using this 

En /7Sical picture, equations (27), (28), amd (29) tollow. Ta 


equation (29), each +t is evaluated at a pressure equal to that 


2/ 





exerted by the gas mixture. Alternately, an expression for 


~ composite would be 


"composite "NON Sr ye 6 


Za aZ ZZ 
where partial pressures are used to evaluate tT for each 
Species. The analogy of the equation for T composite to the 


formula for resistances 1n parallel is apparent. 


B. CHEMICAL NONEQUILIBRIUM FLOW 

In’ order to investigate the effects and mechanics of 
chemical nonequilibrium on the flow in a turbine, the products 
of combustion of a stoichiometric reaction between air and a 
kerosine-type hydrocarbon fuel were considered and a turbine 
expansion calculation performed. The choice of a stoichio- 
metric mixture ratio was indicated by the current trends and 
research goals in the turbojet industry. The internal degrees 
of freedom of all components of the gas were assumed to be in 
equilibrium throughout the expansion. Even though there is a 
known relationship between vibrational modes and chemical re- 
action, this assumption was necessary for analytical BULDOSoen 

The gas was taken aS an inviscid, diffusionless fluid with 
quasi one-dimensional expansion through an ideal turbine. 
The process was considered to be adiabatic. The gas was ex- 
panded first as chemically frozen and then in full chemical 
equilibrium. The stations in the process were as illustrated 
in Figure 8. The concepts used in this section closely follow 


Penner's [13] development for chemically frozen flow.. 
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In chemically frozen flow, the average molecular weight 
obviously remains constant; therefore, the work derived from 


this expansion can be expressed as 


rh 
Work x ee eC ))) (31) 


frozen cal 1 5 


where X. is the molar fraction of species i at station l, n is 
the: number .of.species..present, and H, (T) 1s the total molar 
enthalpy of species 1 at temperature T. 
Given values of PL F T, , and Ty , the total pressure at 
1 ih 2 
Station two was determined from the isentropic relationship 


ys 
yRoL 
jE = Pp —— (32) 


where y* is an appropriately defined average specific heat 


ratio during expansion. Penner [13] gives this ratio as: 


rn 
A ~ eC 
Y* = Sy a a a (33) 
y X, (S°(T, ) - Se (Tt, )) - R &n 2 
i=1 7 a: 2 ee 


The chemical nonequilibrium efficiency, No + was de fined 
as the ratio of the work extracted from chemically frozen gas 
to the work extracted from equilibrium gas expanded to the 


Same total pressure. This efficiency is 


(34) 


Zo 





where 


es Peace ee oe TE) 
bee De ae (36) 
a Dat Ee) tale 2 (37) 
M, = Date ae (38) 
M. - De Pe) (39) 


Once the quantities in the numerator of equation (34) 
were established as functions of known or given quantities, 
attention was focused on the denominator. 

For equilibrium flow, the analysis was restricted to a 
given weight of gas mixture due to the temperature dependence 
of the average molecular weight of the mixture. For purposes 
of comparison it was convenient to choose the initial mole- 
cular weight of the gas mixtures as the characteristic 


weight; thus the ideal or equilibrium flow work was: 


M, 
Work . = HH. - —H (40) 
equi ale M 3 
3 
where 
n 
Hy = z) Xx; (Tt, F Pty) AH + H, (tT, )- H; (t,)| (41) 
i=l 1 1 ue 
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n 
a a jane | aCe eT ) (42) 
3 i=] 1 t. t. | f. vf ty an O 


At this point, all the necessary quantities for evaluation of 


n. were known except T which could be determined from 


Cc e 


3 
isentropic considerations. 

Isentropic flow occurred for the total weight of the gas 
mixture. Since M, was chosen as the fixed weight for analysis, 


the isentropic relationship was: 


S(T oe) (43) 


oes <3 


S (T,P) is the molar entropy of an ideal gas mixture at tem- 
perature T and pressure P and may be expressed as 
= a _ n 
ae ° ——_£ 
S(T,P) = *) X. (T,P)S?(T) R(2n P + a X, (T,P) &nX, (T,P)) 
1=1 i1=1 
(44) 
where S°(T) is the molar entropy at temperature T and one 


atmosphere of pressure. Equations (44) and (44) were combined 


to yield 
n = = ; 
bate, Pe.) sit) —- R &n Pes - R 2&n Xi (Te Pe)| 
= i ¥ x. (2 ,P )|sger )-R &n P 
MeviaT gee: | aS 3 
- R &n a a) (45) 


This equation was used to solve for Ty by iteration. Once 
3 


this has been accomplished, evaluation of equation (34) was 


possible. 
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The mechanics of the chemically reacting flow were then 
investigated. In the combustion chamber of a turbojet, a 
Significant portion of the chemical energy released is used 
to break chemical bonds producing more elementary species. 
This energy 1S unavailable to do work unless it is first re- 
leased through a recombination of the dissociated combustion 
products. The rate of each chemical reaction must be compared 
to a dwell time of the flow at any particular stage in order 
to determine the extent of the unavailability of energy. 


Consider the typical three body reaction 


AtBtmot ape (46) 
*b 

where M is a third body which removes energy during the re- 
combination and rp and rj) are reaction rates in the forward 
(recombination) and backward (dissociation) directions 
respectively. This equation represents two separate reactions 
with the connection that re must equal rh at equilibrium. As 
the temperature and pressure move toward lower values during 
the expansion process, the reaction in equation (46) should 
move in the forward direction. If the reaction rate is not 
fast enough to maintain equilibrium conditions, recombination 
will lag, and the gas driving the turbine will have dissociated 
products whose dissociation energy has not contributed to the 
kinetic energy of the flow. The lag of recombination results 
in a reduction of performance. 


Equilibrium flow calculations involve the simultaneous 


solution of the conservation equations for mass, energy, 


SZ 





momentum, and species, plus the equation of state. This sys- 
tem is simplified for frozen flow because the species composi- 
tion 1S constant. Computer programs for solution of | 
eguilibrium and nonequilibrium flows are presented by Gordon 
and Zeleznik [8] and Zupnik, et al., [23]. 

For relaxing flow in a turbine the problem is complex 
because the nonlinear kinetic differential equations of 
Species concentration and the conservation of atom equations 
are now involved. Solution of the chemical relaxation problem 
is obviously desirable but beyond the scope of the present in- 
vestigation. Penner [12] and Bray [4] have proposed criteria 
to determine when the shift may be made from equilibrium to 
nonequilibrium analysis. Application of these criteria by 
various researchers has produced varied degrees of success. 

Newhall [10] presents twenty-eight elementary reactions 
(Table X) that may be considered significant in the flow of 
combustion products resulting fvoti a mixture of a hydrocarbon 
fuel and air. These equations may be integrated numerically 
with the coupled energy equation to produce concentration- 
time histories for each species during expansion, but this 
was beyond the scope of this thesis. 

The science of chemical kinetics, the study of the rates 
of chemical reactions, 1S complex and difficult, both experi- 
mentally and theoretically. A discussion of both simple col- 
lisional theory and transition state theory is given in Pratt 
[14] and emphasizes the heavy reliance, except in the simplest 


reactions, upon empirical data. 
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The chemical reaction rate may be considered to be pro- 
portional to the total number of collisions between reacting 
chemical species and to the fraction of this total in which 
sufficient energy is available for the reaction to occur. In 
flow systems the rate of reaction is also dependent on the 
mass and energy transport mechanisms. The collision rate is 
dependent upon both concentration and effective collision 
cross section of the molecules involved. 

The most widely utilized approach to the reaction rate 


problem is through the use of the Arrhenius equation 
K-= AG) exp (-E_/RT) (47) 


where A is the pre-exponential factor and E, 1s the experi- 
mental activation energy. The values of both A and E, are 
obtained from spaeeisiwem gal plots of 2n K vs 1/T. Both A and 
E, may be functions of temperature. 

For elementary reactions, such as those that make up the 
complex reactions in tne combustion products, the rate con-. 
Stant, K, from the Arrhenius equation (47) may be employed as 
follows: The rate of change in the concentration of a product 


AB produced by a reaction such as equation (46). 


A+B+Mz?ABt+M (46) 
1s equal to 

el (se) 

Se A) a) (M) (48) 


for an elementary reaction where the parentheses signify 
concentration in moles per cc. 


bed 
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Bach of the equations in Table X had corresponding for- 
ward and backward rate constants from Table XI. Each of 
these rate constants was calculated at the throat of the 
turbine stator model introduced previously (Table XII). The 
mole fractions of the species in the combustion products were 
Peeniculated at the stator throat using Banes, et al., [3] and 
are presented in Tables XIII and XIV. A conversion factor 
that enabled the use of mole fractions instead of molar con- 
centrations in equation (48) was developed and is presented in 
Appendix D. The reaction rates for the equations considered 
were calculated at the stator throat and are tabulated in 
Table XV. The total reaction rate for each specie was deter- 
mined by a summation of the rate of destruction or creation of 
the specie from the applicable equation. If the particular 
Specie appeared on the left-hand side of the reaction equa- 
tion in Table X it was considered as being destroyed at the 
forward rate and produced at the backward. In Table XI, the 
total reaction rate at the stator throat for each specie is 
presented. 

The determination of the presence of chemically frozen 
flow in the turbine required establishment of a characteris-— 
tic reaction time. The total reaction rate, R.y of a specie 
i, had units of moles aS com>, Moen p lyin R; by the 
Specie molecular weight produced units of grams of i per cm? 


per sec, which yields a rate of change of density. 


do. 


oe a 


By taking the gas average density at the throat, Os, and 
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dividing both sides of equation (49) by Ox, a rate of change 
of mass fraction 1s obtained for small mixture density 


changes. 





R, mM; oie dp; : ay. _ 
=— = dt aie 
Ps e) i 
By definition 
ii 


For sufficiently small changes 








dy. sa ama 
eee eae 5 (52) 
dt t, = at 
hetting t, - t, = t, equation (50) may be rewritten as 
i 
4 ee 
R; m5; : ay 7 1. 1 
— =a ie 
oa Si 
Or m m. 
a eee 
2 M, =i 8, 
eR = (53) 
a » oo 
Pi 
Rearranging 
= x x; 
= eal 2 it 
CR om So eae eR (54) 
i i| M M 
2 iL 


Evaluation of equation (54) provided a basis for comparing 


meaction times with dwell times. 
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EEES RESULTS @anD Disetae Le 


Values of the vibrational efficiency defined by equation 
(16) were calculated using stagnation conditions of 2200°K 
and 20 Atmospheres. These values are presented in Table I. 
The data in Table I are plotted in Figure 4 as loss percentage 
against the change of total temperature of the RT gas from 
Stations 1 and 2. See Figure 1 for station numbering. Selec- 
tion of this abscissa illustrates the variance of the vibra- 
tional efficiency with the amount of work extracted. 

The data for different Mach numbers indicates that this 
efficiency as defined is a definite function of turbine exit 
Mach number. It is known that for a given static temperature 
there iS a known amount of energy stored in the vibrational 
mode of internal energy at equilibrium. As the exit Mach num- 
ber increases, the kinetic energy extracts a larger portion of 
the rotational and translational energy depicted in Figure 2. 
Therefore, the static temperature of the RIT gas is lower and 
the amount of equilibrium vibrational energy is lower. The 
total vibrational energy at the turbine exit is the same as 
the equilibrium vibrational energy at the burner outlet. As 
the Mach number increases, the RI gas static temperature de- 
creases and the fraction of the total vibrational energy at 
the exit that corresponds to the RI gas temperature decreases. 
Therefore, the flow experiences a larger departure from equi- 


librium and an attendant increase in the losses. 


bd 
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The apparent minima of curves in Figure 4 as indicated by 
the dashed line is a result of the mechanism of work removal. 
The work was extracted from a RT gas which had both modes 
fully excited in the temperature range considered. This full 
excitation made the work extraction a linear function of tem- 
perature. The equilibrium work in the denominator of equation 
(15) was drawn from all modes of the gas. At the temperatures 
considered, vibrational specific heat is nonlinear and is the 
apparent cause of the minima. 

There were three sources of vibrational relaxation times: 
equation (22), the graphs of Taylor and Bitterman [18], and 
the equations of Anderson, et al., [2]. In Table IV, relaxa- 
tion times for nitrogen-nitrogen collisions are calculated by 
the three methods. The values are for an area change as de- 
fined in equation (26) and stagnation conditions of 2200°K and 
20 atmospheres. The variation in the values of relaxation 
time when calculated by the three methods indicates that the 
data basis for each calculating method is different. The 
graphs of Taylor and Bitterman [18] had insufficient range 
for this analysis. The values obtained through the use of 
equation (22) seem unusually small and therefore were not used. 
The Anderson, et al., [2] data is a later development with 
information for many different collision partners, and its 
use allowed more flexibility of investigation. 

The relaxation times for nitrogen as a function of stator 
area ratio are given in Table V for a stagnation temperature 
of 2200°K and various stagnation pressure ranging from 5 to 


20 atmospheres. Equation (26) was used to calculate area ratio 
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as a function of distance along the midspan streamline. These 
data are plotted in Figure 5 with the area ratio as the 
abscissa. The area ratio values are of significance only in 
that they define data collection points. The pressure de- 
pendence of the relaxation is aptly demonstrated in this 
figure. The lower pressure results in a lower collision fre- 
quency and thereby slows the relaxation. The curve of relaxa- 
tion time for nitrogen-water collisions emphasizes the 
catalytic effect of water on the relaxation process. The 
presence of water in the combustion products of a hydrocarbon 
and air mixture helps to delay the onset of vibrational 
freezing. 

The dwell time as plotted in Figure 5 represents the re- 
ciprocal of the local velocity in units of seconds per inch 
multiplied by a factor of ten. The order of magnitude of the 
dwell time was raised because exact numerical solution of 
equation (21) indicated vibrational freezing is Significant 
at A/A* = 1.14. The factor of ten multiplication results in 
the dwell time and relaxation time curve intersection agree- 
ing with that solution. The multiplication factor makes the 
dwell time longer and more pessimistic. The intersection of 
the dwell time curve with the relaxation time curves represents 
a position downstream of which the flow may be considered to 
be vibrationally frozen. For the pure nitrogen, this point 
occurs at a low area ratio indicating that the flow is prob- 
ably fully frozen for practical calculations. Once the flow 
has frozen, Bray [4] indicates that it will remain frozen or 


at least in some form of nonequilibrium. The fact that 
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freezing is indicated early in the stator section of the 
turbine implies that the flow should be considered as vibra- 
tionally frozen throughout the turbine performance 
calculations. | 

Equation (21) was evaluated for various conditions and the 
data plotted in Figures 6 and 7. Both of these figures dra- 
Matically demonstrate the departure from vibrational equilib- 
rium that occurs at a low area ratio. The effects of other 
collision partners are shown in Figure 7 using data from 
Tables VII and VIII. The upward movement of the curve in 
Figure 7 when the midspan streamline length 1s reduced is a 
result of a lower gas transit time. 

Values of the chemical nonequilibrium efficiency were com- 
puted from equation (34) for stagnation conditions of 2200°K 
and 20 atmospheres of pressure at the burner outlet. These 
values are presented in Table IX. Data from Table IX were 
plotted as loss percentage in Figure 9 with the change of 
total temperature from burner outlet to turbine exit as the 
abscissa. Here again the total temperature change is an indi- 
cator of the amount of work extracted. The losses due to 
chemical freezing in this stoichiometric-mixture expansion are 
rather low. For the temperature range considered, the amount 
of nitrogen dissociation is negligible. This is significant 
because nitrogen comprises nearly three-quarters of the gas 
mixture. For the stoichiometric mixture there is very little 
oxygen remaining to dissociate after the reaction. Consider- 


ing the size of the losses, the tendencies of the curve are 
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en oce to question. The precision of the thermochemical 
data utilized for the calculation would have significant ef- 
fect on the results. 

The characteristic reaction time for each specie was 
calculated from equation (54) using the data for the stator 
throat and stator exit at total temperature of 2200°K and 
total pressures of 15 and 5 atmospheres. These times are in 
Table XVIII. Only the reactions considered to be the fastest 
were utilized in the analysis. The influence of slower reac- 
tions upon the fast reactions was not analyzed. Each of the 
times was analyzed and a classification was performed. The 
gas transit time from stator throat to exit was compared to 
each specie reaction time, and the specie was placed in an 
appropriate section of Table XIX. 

A negative reaction time for a Specie indicated that the 
predominate or fast reactions involving that Specie were away 
from equilibrium for the particular conditions of temperature, 
pressure and concentrations. The real life case, including 
reactions not considered here, would, for long time periods, 
cause the reaction times to become positive. The negative 
reaction time species were classified as frozen. A long or 
Short positive reaction time compared to the transit time 
classified the specie as frozen or equilibrium respectively. 
Specie reaction times of the same order of magnitude as the 
transit times were classified as unknown. The approximation 
used to calculate the reaction times was not sufficiently 


accurate to make a definite statement as to the freezing of 
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the particular specie. Table XIX definitely indicates the 
presence of chemically near-frozen flow. 

Several aspects of the current state of the art in turbo- 
jet development tend to reinforce the gravity of the losses 
calculated in this paper. As the aircraft operate at higher 
and higher altitudes, the total pressure at the burner outlet 
will decrease. This will insure the presence of vibrational 
and chemical nonequilibrium as illustrated in Figure 5, 6, 7, 
and Table XIX. The use of hydrogen fuel would produce large 
amounts of water thereby delaying the vibrational freezing 
effects. The magnitude of the burner outlet temperature will 
vary as the ability of the turbine blade cooling method to 
handle high inlet temperature, and this variation will effect 


the degree of nonequilibrium present. 
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section 
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The 


IV. CONCLUSIONS 


vibrational freezing of the flow in the turbine 

of an advanced turbojet using a stoichiometric mix- 
kerosine and air can result in losses of up to Seven 
when compared to equilibrium flow. 


Vibrational freezing of the flow will be virtually 


complete a short distance after the throat of the first 


stator assembly. 


The 


losses due to chemical freezing are insignificant 


compared to the vibrational freezing losses. 
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APPENDIX A 


TABLE I 


VIBRATIONAL EFFICIENCY WITH TOTAL TEMPERATURE OF 2200°K 


AND TOTAL PRESSURE OF 20 ATM 
t ° fig: 1 — 2 
mt, °K Mt = .1 Ms = .3 es 
1500 96.37 95.21 32.86 
1550 96.85 95.40 PGs 
1600 97.00 95.52 92.49 
1650 97.27 95.95 aes 
1700 97.71 95.89 92.02 
1750 97.90 95.80 iL 62 
1800 97.98 95.80 91.03 
1850 97.81 95.43 90.51 
1900 98.14 95.14 89.46 
TABLE II 
VALUES OF CONSTANTS FOR USE IN EQUATION (22) 
FROM VINCENTI AND KRUGER [20] 
; Heat Bath c e Approximate Range 
Species | wolecules|Atm-microsec| a Of Kk 
0, 0, ped2 2g) ongS eo” Sevan e700 
ey Ar B85 x) 10) 95a ee) s0cm= 4300 
N, N, 7.12 x 107) 1-91 10°) 9800 ="6000 
NO NO me6 x 10) 1? x 10" «ASO Om= 8000 
NO NO pelo x 10m =e? x 1028 1S0Cmemtoan 
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TABGE iil 


STATOR AREA CHANGE WITH ASSOCIATED IDEAL GAS DATA 


FROM WANG, ET AL. [21] 
y = 1.26 
Velocity, Dwell Time, 
* 
A/A a wy eee M/sec sec/in 
lL. weleeal) - 885 mes eal 854 DoT 

1.14 eas 219 ol 1154 22 
me2/ 1.59 57526 ~2523 252 20% 3 
P37 1.69 e722 wale, T3200 19.4 
1.46 7 a7 GG ail 2 as 1354 18.8 
Pol inet s 702 nS 1376 £325 
wo ios -696 pelegZ 1389 Ee. 3 

TABLE IV 


COMPARISON BETWEEN VARIOUS METHODS OF OBTAINING STATOR 


RELAXATION TIMES WITH TOTAL TEMPERATURE = 2200°K 
AND TOTAL PRESSURE = 15 ATM 
A/as T P Tegn. (22) ‘Graphs ‘Eqn. (23) 
1 1947 11.06 re 3 635.5 Gilaul! 
| 1740 6.3 38. 7 172.4 207-5 
27 1656 5.05 50.6 238 353 
is 3 / 1604 4.34 65.8 Ke 498 
1.48 1563 33°65 Silas kk - 658 
ee. 1544 3.60 eile [8 Ke 759 
ie 55 1531 Sale! 98 K* 837 
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** Unable to extrapolate the graph further. 





A/A* 


TABLE V 


STATOR RELAXATION TIMES IN MICRO-SECONDS 
OF NITROGEN MOLECULES UNDER VARIED TOTAL PRESSURES 
WITH TOTAL TEMPERATURE 


20 


61.1 
20725 
353 
498 
658 
7159 


837 


15 


82.9 


277 
471 
663 
876 


1011 


"1116 


TABLE VI 


2200°K USING EQUATION 


Pressure, 


15 


*k* 


NOs: 


2. 
28 
35 
Bg 
42 


44 


10 


124 
415 
706 
996 
1316 
1577 


1673 


*k Collision Partner is H,O in this column vice N 


as in the other “columns. 


(23) 


COMPOSITE RELAXATION TIMES CALCULATED USING EQUATION (30a) 
WITH TOTAL PRESSURE = 


A/A* 
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jg On 


Z 


10.8 
Za 
Z2Tsd 
33.4 
BiOar! 
42 


44.3 


15 ATM 


2 


Composite 


44.2 
106.4 
2 2) 
190.6 
228 
250.5 
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RATIO OF ACTUAL VIBRATIONAL ENERGY 
TO EQUILIBRIUM VIBRATIONAL ENERGY AT THE THROAT 


TABLE VII 


UNDER VARIOUS TOTAL PRESSURES AT TOTAL TEMPERATURE = 


Total 
Pressure, 


20 


cS 


ATM 


E* 
V 


E* 





A/A* 
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Total 
Pressure 


10 


A/A* 


2200mi 








TABLE VIII 


RATIO OF ACTUAL VIBRATIONAL ENERGY TO EQUILIBRIUM VITRATIONAL 
ENERGY AT THE THROAT WITH TOTAL PRESSURE = 15 ATM AND TOTAL 
TEMPERATURE = 2200°K FOR VARIOUS MIDSPAN STREAMLINE LENGTHS 














EX E., EE. EY 
Ex A/A* Ex Rx loka 
Yth “th Yth Yth 
ii Wr Tae lr 
i i 1 : 1 iL 
-8388 1.14 -9938 72979 -9946 
m7 733 iee7 91832 ~9942 -9846 
m7 3 1 Leoe oT oak ~9914 29 fois 
- 7706 1.46 -9687 Serial - 96 80 
-6856 le 5 29633 ist -9609 
BO7 95 o> -9584 ~9854 ~9544 
t 3: Midspan streamline length is 6"; collision Partner 
N.- 
tt : Midspan streamline length is 2"; collision Partner 
No- 
ttt : Midspan streamline length is 2"; collision Partners, 
Na7eco. ane 1-0. 
2 2 zZ 
TABLE IX 


CHEMICAL EFFICIENCY CALCULATED USING EQUATION (34) 
WITH TOTAL TEMPERATURE = 2200°K 


Turbine oxi 





. Total Temperature @ 
1500 513) 41 
1600 99.00 
1700 2) oe) lls: 
1800 99515 
1900 99 0 
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ELEMENTARY CHEMICAL REACTIONS FROM NEWHALL [10] 


NO 


NO 


OH 


OH 
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TABLE X 


ty ty t+ ty ty hy t+ i ay ty ty¥ ty¥ ty¥ 


ty 


OH 





TABLE XI 


SELECTED RATE CONSTANTS 


FOR 


Forward Rate Constant 


THE ELEMENTARY CHEMICAL REACTIONS FROM NEWHALL [10] 


Backward RAte 


Reaction Concune 
1 5.4 x 107’ exp(-123,600/RT) 1.5 x i0%° 
2 3.1x10*> exp (-110,000/RT) 7On, Mie 7 
10 2 
3 5 C150 000) exp(-150,000/RT) 9.0 x 1074 
T (RT) 
4.2x107* 224,900 14 
4 - T (————) exp (-224,900/RT) 6 See 10 
, 2 RT 
5.4x107+ iG 
5 = exp (-74,000/RT) 2.05% 10 
6 6x10+> ane wet exp(-11,800/RT) 1. On x Lae" 
IL 12 
vi 1.4x10° "exp (-6,000/RT) 3.3x10°- "exp (-8,000/RT) 
13 14 
8 5.5x10° ~exp(-1,000/RT) 7.2xX10”~ exp (-16,900/RT) 
ILS 14 
9 6.2x10°~exp(-6,000/RT) 3.2*1075 sexp (—-2 Pane Re 
12 13 
10 7.7x10 “exp (-1,000/RT) 8.3x10° “exp (-18,100/RT) 
ip 14 
‘hal 7.1x10° “exp (-7, 700/RT) A. fX1O = sexo (-2 772507 Re 
12 3.2x10" exp(-39,100/RT) 13.3x10%exp (-7,080/RT) 
13 Tescesia = 7x10+%exp (-75,500/RT) 
14 1.8x107° D4 2 exp (-47,000/RT) 5.ex10l9 pl/2 
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Reaction 


13 


14 


TABLE XII 


RATE CONSTANTS AT 1947°K 


Forward Constant 


RIG FO) a ore 


1.3896 x 10° 


SLSR sls 


Boe se 10) dae 


1.36733 x 102° 


1.037 x 10° 


2.9688 x 107+ 


4.247 x 102° 


jest coe 


5.9461 x 10° 


9.702 x 1ott 


2.542 x 10° 


1.55 x 1023 


4.205 x 10° 


oye 


Backward Constant 


1.5 x 102° 


3.5933 x 1017 


9.0 x 102° 


6.1 x lot? 


2.0 x 102° 


i ps 0 a’ 1017 


4.17295 x 107+ 


9.1234 x 10° 


1.3693 x o> 


‘lees <2 eo 
4.1025 x 10 
4.1536 x 10 
2.3418 x 10 


2.559 x 102° 





TABLE XIII 


EQUILIBRIUM MOLE FRACTIONS WITH TOTAL TEMPERATURE = 2200°K 
AND TOTAL PRESSURE = 15 ATM FROM BANES, ET AL., [3] 








Throat Jp. gal os 
T 1947°K 1531°K 
P Se 6) shell 2 eo Em 
CO, 1.2919x1077 1.3073x1l0 > 
H,0 1.2992x107+ 1.3048x10~+ 
0, 6.7092x10 4 5.3073x10 > 
H, 3.4679x10 * 3.6215x10 > 
N, 7.2896x10-* 7.2981x10 > 
N aoa 1x107)4 
NO 3.8768x1l0* 2.4451x10~> 
OH 2.5755x10~- 1.2768x10~> 
CO 1. 4098x1077 1.0068x10* 
NO, 1.1142x107/ 2.8080x10~” 
H 7.9334x10~° 1.2395x1o07’ 
O 4.4208x10 ° 3.8529x10~° 
M 28.876 28.903 
5 1.499x10~° 5.938x10- 
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TABLE XIV 


EQUILIBRIUM MOLE FRACTIONS WITH TOTAL TEMPERATURE 
AND TOTAL PRESSURE = 5 ATM FROM BANES, ET AL., 


=e a 


so} 


Throat 


1947°K 
2.7655 atm 
1.2845x10_ 
1.2963x10> 
1.0156x10 
5.0574x10- 
7.2857x10- 
10 tae 
4.7613x10_ 
3.8269x10> 
2.1936x10 
9.4140x10 
1.740ex 10" 
9.8841x10 
28.863 


4.997x10 4 
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al 


a 


3 


4 


iL 


4 


4 


3 


8 


3) 


6 


Exit 


1531°K 
0.86 atm 
1.3068x10 
1.3045x10- 
7.6494x10- 
5.0628x10_ 
7.2979x10_ 
i197 
2.9394x10 
e210 
1.7399x10 
2.4052x10> 
2.4370x10> 
7.6885x10 
28.902 


1.98x10~7 


—_— 
—_ 


TL 


i. 


5 


5 


u 


5 


5 


4 


9 


7 


8 


2200 K 
[3] 





_—— 


Total Pressure, 


aol: 


Forward 
1.839x10 ° 


9.48x10 1° 


2.636x10 2! 


8.748x10 2? 


2.997x10~° 


1.369x10 2° 


1.637x10 ° 


1 SOG 


3.169x10 > 


1.064x10~° 


1 OUEST 


oie © 


1.620x10 11 


2.952x107? 


TABLE XV 


iS 


Backward 


Aes 8 


23428 210 


4.066x10 ~° 


1.02x10~2/ 


3.503x10 ° 


1.997x107 29 


1.727x10~° 


Te elie 


3.809x10 > 


LSBs a ae 


1.1348x10 > 


7.52x10 1 


2.04x107? 


3.40x10~° 
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REACTIONS RATES OF ELEMENTARY EQUATIONS AT 


Forward 
2.037x10 


1.534x10 1° 


3.593x10 28 


9.703x10 7° 


2 esi © 


2.30x10 21 


5 Sosg ” 


aie emia 


7.62x10~° 


2.61x10~4 


Ace 


3.58x10 2° 


2.21x10 1 


6.09x10~19 


1A Teak 


5 


Backward 
3.775x107 / 


4.114x19~2° 


3. 36x10 1° 


2.305x10 7% 


3.56x10 7! 


3.6910 


eal 


Rese 


9 26210 


2.28x10 > 


2.75x10° 


1.26x10 > 


5.05x10 29 


7.14x107 19 





TABLE XVI 


EQUATIONS FOR TOTAL REACTION RATE OF EACH SPECIE 


Selo), aie Salis 

HO ssp reakiek s0\sem—- Cl tae oan ils 

*o,, “6p ‘of (8f “8b “12£ "12b°“l4b “14£ 

“aL fob (2f “TE 7b Oe 

“N5 Yap *4£°°13£ 713b 

Ryo Mab (3£ “Sf “Sb ob -12f “1b "136 “aape Lar 

Ron “1£*1b**7b-*7£"* gb“ 8f* "9b "9f°**1ob **10f"11b *11£ 
Roo eel iene sate 

Ryo (6b SE o14t Ab 


1 ine ee in Oe ee ace ae 


-R., +2r-,.--2 


Be 3c ob Nee Se Cre Ch 7fe 7b cee oe Oo 


[pr lo ee eee ea ee Sala 
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TABLE XVII 


TOTAL REACTION RATE, R, AT THE THROAT 
WITH STATIC TEMPERATURE = 1947°K 


Total Pressure 


Species 
15 5 
-4 as 
CO, Seco = yee 
-4 -3 
H,0 = Gaal ~ 2.185x10 
=F -7 
07 - 4.,98x10 - 1.695x10 
H, 6.4x10 - 1.63x10 4 
N, sain ~ 5,03x10 1° 
~7 -8 
NO - 5.047x10 - 7.445x10 
, <8 <8 
OH 1.019x10 4.24x10 
-4 Ss 
co iL AUS EEE 3.13x10 
~7 -g 
NO, 5.056x10 7.46x10 
H _ 9 Bosca = iL Selo 
A -3 
0 = 1 Sieh ~ 2.022x10 
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TABLE XVIII 


CHARACTERISTIC REACTION TIMES IN THE THROAT AREA 
IN ORDER TO REACH EXIT EQUILIBRIUM CONDITIONS 
FOR TOTAL TEMPERATURE = 2200°K 


Total Pressure 


IS 5 
Species Time 
CO, fs (R. ee = + 1.14 
ae 3 
H.0 raf ane eee 16 
0, = Gadds 10 oe ~ 959x107 
ae 5) 
H + 2.52x10 + 4,.84x10 
N, + 4.32 +e eineeto 
NO ~ 3.74x10~7 ~ 1.04107 
<5 = 
OH + 1.25x10 + 1.49x10 
4 1 
CO + 6.26x10 + 2.00x10 
es <1 
NO. Teleaioe: ia fey Goan 
H = eal 5 A seo 
0 Evy oe = Ge sonia 
leases lee Be eclome 
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TABLE XIX 


FREEZING CLASSIFICATION OF CHEMICAL SPECIES 


Pressure 


15 None 


None 
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APPENDIX 8B 


Stator Exit Guide Vanes 
COMBUSTION | | || | a 
tt 
ae 


FIGURE | TURBINE SCHEMATIC AND FLOW STATION 
[IDENTIFICATION FOR THE VIBRATIONAL 


FREEZING ANALYSIS 


oo 





Equilibrium 







Total 


/ 


Work Extracted 


2 3 AXIAL 


DISTANCE 
Frozen Pn en a 


FIGURE 2 ENERGY DISTRIBUTION IN THE TURBINE 
FOR THE VIBRATIONAL FREEZING ANALYSIS 
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FIGURE 3 


STATOR PASSAGE 


Midspan Streamline 





lO 
a © 
O EOCwS OF 
_j S MINIMA 
6 / 

iu 8 - 
= © S us @ 
Zz 1 
<[ y M, = 5 
= / 
Oe / 
ui 6 / 
nt / 
a. 4 

> vi M= 3 
ul) S 4 O © 
O) U S ® 
a+ 
= wa 
= ae 6 M. = 
Lil ee (-) 2. 
eae 
LJ Os @ 

300 400 500 600 (00 


CHANGE IN TOTAL TEMPERATURE ACROSS THE 
TURBINE, Ty, - Ty,, °K 


FIGURE 4 PERCENTAGE PERFORMANCE LOST DUE 
TO VIBRATIONAL FREEZING 
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5 AT 


Na Ne 10 AIM 


lOO0O SS 15 ATW 
\ si 


= 


Dwell Time 
lOO?F 
W) 
ra) 
S 
O I5 ATM Neo-H OQ 
lJ 
7p) 
© 
a 
© 
= 
A lO 
= 
b= 


| 1.14 |.27 [37 1.46 [5] 
AREA RATIO, A/A* 


FIGURE 5 LOCAL VIBRATIONAL RELAXATION TIMES FOR 
: THE STATOR PASSAGE IN FIGURE 3 
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98 
O 
96 © ID ATN © 
20ATM C 
94 
Ey 
aa 
=> th 
No-N, COLLISIONS 
MIDSPAN STREAMLINE LENGTH 6° 
90 
88 
Equilibrium 
86 
| 114 |.27 1.37 1.46 15I 


AREA RATIO, A/ A* 


FIGURE 6 RATIO OF LOCAL VIBRATIONAL ENERGY TO 
EQUILIBRIUM THROAT VIBRATIONAL ENERGY 
AS A FUNCTION OF AREA RATIO FOR 
PURE NITROGEN 
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98 ——— 


o 


SYMBOL MIDSPAN STREAMLINE COLLISION 


LENGTH PARTNERS 
94 
A 6" No-~No 
Ey 
Es 0) 2" No - No 
th 
92 
a 2 No -No 
No —H»O 
90 
«—— EQUILIBRIUM 
88 
86 
1.14 1.27 137 146 |.5| 


AREA RATIO, A/A* 


FIGURE 7 INFLUENCE OF STATOR SIZE AND COLLISION 
PARTNER ON THE RATIO OF LOCAL VIBRATIONAL 
ENERGY TO EQUILIBRIUM THROAT VIBRATIONAL 
ENERGY AS A FUNCTION OF AREA RATIO 
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4 | 





COMPLETE EQUILIBRIUM 





ENTROPY 


CHEMICALLY FROZEN 





ENTROPY 


FIGURE 8 EXPANSION PROCESSES USED _ IN 
ANALYSIS OF CHEMICAL FREEZING 
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PERCENTAGE PERFORMANCE LOST 


300 


CHANGE IN TOTAL TEMPERATURE ACROSS THE TURBINE 
-T, , °K 


OW tee 





400 900 600 00 


PERCENTAGE PERFORMANCE LOST AS A 
RESULT OF CHEMICAL FREEZING 
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APPENDIX C 


An approximate solution to equation (21) may be obtained 
using the finite difference method. Equation (21) is put in 


the form as in Hurle, et al., [9]. 


[l - (Ax/2A_)] 


BCP nen > Tr ax7ax_ pT Pv) n 
: Ax/2 BN) 5) r ECT) im 
aa 2 ee ve 


where n and n+ 1 refer to two adjacent locations along the 
nozzle axis Separated by distance Ax, and A and T are evalu- 
ated for equilibrium flow. A is the characteristic relaxation 


length and »} = Ut, where U is the local gas velocity and tis 


relaxation time. 
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APPENDIX D 


This is the calculation of the conversion constant for 
use in the rate equations so that mole fractions may be 
utilized directly. 

The units of the rate constants in Table XI vary with 
the type of elementary equation; bi-molecular or tri- 
molecular. The normal units for concentrations are moles 


3and the units of X are moles of specie 1 per mole of 


- cm 
mixture. From the ideal gas law 
POV =en Rk T 


Multiplying by moles of specie 1 Over moles of specie 1 





PV X. = n. RT 
ne 7 
or 
Eee 
V R T 
in units of moles - Be the desired units. 


Using the gas constant equal to 


oe atm 
mole °K 


R = 82.06 
a conversion factor for use at a specified temperature and 
mixture pressure can be obtained. The mole fractions are each 


multiplied by this factor and are then ready for use in the 


Specie rate equations. 
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